Emergence of continual directed flow in Hamiltonian systems 
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We propose a minimal model for the emergence of a directed flow in autonomous Hamiltonian 
systems. It is shown that internal breaking of the spatio-temporal symmetries, via localised initial 
conditions, that are unbiased with respect to the transporting degree of freedom, and transient 
chaos conspire to form the physical mechanism for the occurrence of a current. Most importantly, 
after passage through the transient chaos, trajectories perform solely regular transporting motion 
so that the resulting current is of continual ballistic nature. This has to be distinguished from the 
features of transport reported previously for driven Hamiltonian systems with mixed phase space 
where transport is determined by intermittent behaviour exhibiting power-law decay statistics of 
the duration of regular ballistic periods. 

PACS numbers: 05.60.Cd, 05.45.Ac, 05.60.-k, 05.45.Pq 



Nonlinear transport processes of particles evolving in a 
spatially periodic potential have attracted considerable 
interest [1| (for a recent review see In most of the 
studies, the emergence of particle current is triggered by 
an external time-dependent field with zero mean (which 
can be of stochastic nature) or is provided by a determin- 
istic periodic force. Recently, the Hamiltonian case has 
become the subject of intensive studies due to its rele- 
vance to the motion of cold atoms in an optical potential 
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The necessary conditions for rectification of the cur- 
rent, based on symmetry investigations of the external 
field and the underlying static potential, have been pre- 
sented in [1] and To be precise, all symmetries that, 
to each trajectory, generate a counterpart moving in the 
opposite direction, need to be broken. This is achievable 
by imposing a time-dependent external force that is peri- 
odic but not symmetric under time reversal Fur- 
thermore, the phase space has to possess a mixing prop- 
erty, with coexisting regular and chaotic dynamics [6|. 
In extended chaotic systems a nonzero current can be 
obtained as the time-averaged velocity of an ensemble 
of trajectories in the chaotic component of phase space, 
and the chaotic transport proceeds ballistically and di- 
rectedly @, [|. 

Extensions to studies of autonomous Hamiltonian sys- 
tems of one-dimensional billiard chains have followed |9j , 
[ioj . The necessity of creating chaos requires at least two 
degrees of freedom. As an example of such a system, a 
classical magnetic billiard for particles carrying an elec- 
tric charge has been studied in Q. In order to break 
the time-reversal invariance, an external static magnetic 
field, penetrating the plane of motion perpendicularly, 
has been applied. In addition, achieving directed trans- 
port demands breaking the remaining spatial symmetry, 
which can be achieved, e.g., by properly placed asymmet- 
ric obstacles inside the billiard @, [l(|. Uni-directional 
motion in a serpent billiard chain has been reported 
in [Hi. 

The aim of the current work is to demonstrate that 
chaotic directed transport as achieved in systems with 



a mixed phase space is not the only option to obtain 
directed transport in Hamiltonian systems. In driven 
Hamiltonian systems the directed transport necessitates 
a mixed phase space for which chaotic trajectories can 
stick to the boundaries of regular regions inducing long 
periods of nearly regular motion so that finite asymptotic 
currents can be observed On the other hand, 

sustained transport seems to be impeded by the inter- 
mittent behaviour exhibiting power-law decay statistics 
of the duration of periods of regular motion. Therefore, 
for the sake of stability and reliability of a transporting 
regime it is desirable to find ways to accomplish contin- 
ual regular transport. As we show, when dealing with 
nonintegrable systems, it is advantageous when chaos is 
only of transient nature and serves to guide trajectories 
onto regular transporting motion. Furthermore, we also 
show that, in contrast to the studies quoted above, a di- 
rected flow can arise in autonomous Hamiltonian systems 
even without the application of a time reversibility sym- 
metry breaking external field: specifically, we show that 
while the system as a whole is time-reversible, physically- 
relevant sets of localised initial conditions lead to current 
despite being unbiased with respect to the transporting 
degree of freedom. 

We consider a Hamiltonian system with n degrees of 
freedom of the form 



H{p, q) = ^p 2 



(1) 



with (p, q) £ R 2n and U (q) as a potential function. It is 
assumed that the system possesses an open component, 
by which we mean that constant energy surfaces may be 
unbounded in the coordinate(s), allowing for transient 
chaos. 

In the following, we discuss the spatio-temporal sym- 
metry properties of the Hamiltonian equations p — 
—dH/dq and q — dH/dp. Firstly, the system of equa- 
tions exhibits time reversible invariance. Solutions are 
of the form X(t) = (p(t),q(t)). Applying the time- 
reversal operator yields r(p(t),q(t)) = (—p(—t),q(—t)) 
and, hence, if X is a solution, then so is tX. As for 
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the implication of time-reversibility symmetry with re- 
gard to the net flow, let a solution, starting from some 
initial condition X(0), be evolved in time up to a finite 
observation time T at which the forward trajectory ar- 
rives at the point X(T) on the constant energy surface. 
Subsequently, under application of the time-reversal op- 
eration, the signs of the momenta at this point, X(T), 
are reversed, and letting the solution evolve once again, 
with tX(T) as the initial condition, the corresponding 
backward trajectory traces back the path of the forward 
trajectory in coordinate space. For a microcanonical en- 
semble, the initial conditions X(0) and tX(T) are equally 
selected points from the constant energy surface. Thus, 
for systems with time-reversibility symmetry and uni- 
formly distributed initial conditions populating the whole 
energy surface there is no preferred direction of the flow 
preventing the emergence of a current. 

However, energy surfaces that are unbounded (along 
the coordinates of the open component) may not, in prac- 
tise, be completely populated with a finite set of initial 
conditions. In this sense, any finite set of initial con- 
ditions can be regarded as being localised in the coor- 
dinates of the spatially-open system. Localised initial 
coordinates are frequently used in applications such as 
for the problem of a particle flow emerging when the 
particles are initially trapped in a well of a spatially 
infinitely extended, multiple- well potential (see below). 
Therefore, for practical purposes, it is supposed that the 
initial conditions for coordinates are localised in the do- 
main — q.jj < qj(0) < qj.r and 1 < j < n. We define 
a trajectory of the spatially-open system as transporting 
(with respect to a specified observation time, T) if: (i) at 
least one of its coordinates, qj(t), escapes from the do- 
main of the localised initial conditions at some instant of 
time < t escape -C T and, (ii) performs subsequently di- 
rected net motion, that is (jPjit)) ^ for T >t > t escape , 
where ( ■ } denotes the time average. This leads the tra- 
jectory away from the domain of localised initial condi- 
tions so that at the end of the observation time, at least 
one of the terminal coordinates of a transporting forward 
trajectory lies outside the domain of the localised initial 
values; qj{T) < —qj t i or qj(T) > qj tr . Consequently, the 
initial condition of the corresponding backward trajec- 
tory (which would compensate the contribution of the 
forward trajectory to the net current) is not contained 
in the set of localised initial conditions. It needs to be 
stressed that this alone does not imply the emergence of 
a current in the system for such sets of initial conditions. 

Furthermore, although the equations of motion are 
time reversal symmetric, not all of their solutions nec- 
essarily obey this symmetry. Regarding the selection of 
an initial condition, X(0), it holds that unless a trajec- 
tory is self-reversed, i.e., for the selected initial condi- 
tion and its time-reversed counterpart, tX(0), the corre- 
sponding motion coincides in phase space, i.e., X = tX, 
time reversibility is broken. In fact, selecting an initial 
condition, X(0), may break the time-reversibility sym- 
metry, which is the case when the trajectory X and its 



time reversal counterpart tX are distinct. Examples 
of time-reversibility breaking are provided by the (un- 
bounded) rotating trajectories of a pendulum whereas 
the (bounded) librating trajectories are self-reversed and 
thus time reversible symmetric. For potential systems of 
the form of Eq. |T]) , time- reversibility is manifested in co- 
ordinate space in the spatial symmetry features induced 
by reflections on the time-reversal symmetry hypersur- 
faces on which the time evolution of the trajectories sat- 
isfies p(— t) = — p(t) and q{— t) — q{t). They are repre- 
sented by (n — l)-dimensional manifolds in the coordinate 
space and are determined by the condition p k = 0. This 
gives 

dU 

S k : — = F k (q)=0, l<k<n. (2) 
oqk 

(For the sake of illustration, we consider the case where 
each condition in ^ has a single solution and the re- 
sulting symmetry manifolds intersect in a single point 
representing an equilibrium of the system.) Let us con- 
sider reflections of a trajectory, projected onto coordinate 
space, on the symmetry manifolds S k induced by the cor- 
responding operators Rk- Self-reversed trajectories are 
left invariant upon reflections on the symmetry mani- 
folds, rifc=i Rk(q) = Ii since they are mapped pointwise 
onto themselves on equipotentials U(q) = U(Rk{q)) such 
that the sign on the r.h.s. in the equations of motion 
for p k is reversed, F k {q) -> -F k (R k (q)), 1 < k < n. 
Observe that upon reflecting on all symmetry manifolds 
the relation U(q) = U(Y\™ =1 R k (q)), 1 < m < n is 
left invariant under permutations of the reflection op- 
erators. In fact, time-reversing symmetry is in coordi- 
nate space tantamount to invariance with respect to re- 
flections on the symmetry manifolds. In more detail, 
any self-reversed trajectory, projected onto coordinate 
space, repeatedly crosses every symmetry manifold S k 
each time with an opposite sign of the corresponding 
force — F k (q) < oo. Moreover, successive crossings of 
a single symmetry manifold, S k , occur in alternating di- 
rections. Thus there must be turning points for the tra- 
jectory implying bounded motion and no directed flow 
can arise. Notice that no assumptions with regard to 
the spatial symmetries of the trajectory are needed. In 
contrast, as transporting (unbounded) trajectories are 
not invariant with respect to reflections on the symme- 
try manifolds, preservation of time-reversing symmetry is 
not possible for such a single trajectory. A transporting 
trajectory may escape without having crossed a symme- 
try manifold at all. However, if it does cross then after 
all such crossings of a symmetry manifold, the escaping 
trajectory promotes directed transport. Nevertheless, re- 
flections on the symmetry manifolds, mapping a trans- 
porting trajectory onto another transporting one, can in- 
duce spatial symmetries such that these two trajectories 
mutually compensate each others contribution to the net 
flow. Let the point qo in coordinate space be an ini- 
tial condition associated with a transporting trajectory. 
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Reflecting in coordinate space on the symmetry mani- 
folds Sk transforms an original point, qq, into its image 
point, qi, according to Rk<lo = Qi,k- While the value 
of the potential energy is maintained, U(qo) = U(qi^ k ), 
the sign on the r.h.s. in the equations of motion for 
Pk is reversed; F k {R k qo) -> -F k {qi. k ). However, the 
magnitude of the gradients, F k — dlf/dq k , is not nec- 
essarily maintained. Reflection on all of the symmetry 
manifolds yields nSb=i R-kqo = Qii reversing the sign in 
all of the r.h.s. of the equations of motion for the evo- 
lution of the momenta F k {qj) — > —F k (qo),l < k < n. 
With the time evolution of a coordinate expressed as 

q(t) = q(0) + f* dt'{p(0) + /*' dt"[-F(q(t"))]},we con- 
elude that, for the pair of trajectories emanating from 
qo and qi, symmetry (zero net flow) results if (pi,qi) — 
(~Po, -Qo) so that F k (q r ) = -F k (q Q ), 1 < k < n. This 
is the case when the potential is even in the coordinates, 
that is U{q) = U(—q). Then there exist pairs of current- 
annihilating counterpropagating trajectories, X(t), start- 
ing from X(0), and — X(t), starting from — A(0), re- 
spectively. In other words, reversion symmetry under 
reflections on the symmetry manifolds is needed for zero 
net flow which, together with invariance with respect to 
changes of the sign of the momenta, amounts to parity- 
symmetry of the system H(p,q) = H(—p,—q). Con- 
versely, violation of reversion symmetry with respect to 
at least one of the coordinates q k establishes a prerequi- 
site for the occurrence of directed flow. 

For an illustration, we consider the conservative 
and deterministic dynamics of a particle whose coordi- 
nate q evolves in a one-dimensional periodic, spatially- 
symmetric washboard potential of unit period which is 
given by 



system of coupled equations of motion is given by 



U(q) = U(q + l) = 



1 — cos(27rq) 
2^ 



(3) 



The particle is assumed to interact with a local oscillator 
of amplitude Q evolving in a harmonic potential being 
subjected to a tilt force of strength F 



1 



V(Q) = -u 2 Q 2 + FQ 



(4) 



Coupling between the particle and the harmonic oscilla- 
tor arises from the interaction potential 



W(q,Q) = D 



1 - 



1 



cosh(<7 — Q) 



(5) 



that depends on the relative distance \q— Q\. The param- 
eter D regulates the strength of the coupling between the 
particle and the harmonic oscillator degree of freedom. 
Most importantly, the interaction is of local character, 
having negligible influence on the motion outside the so 
called interaction region, in which when the relative dis- 
tance \q — Q\ is sufficiently small so that a significantly 
large gradient of the interaction potential results. The 



q = — sin(27r<7) — D 



tanh(<7 — Q) 



-Q-F + D 



cosh(g — Q) 
tanh(<7 — Q) 
cosh(<7 — Q) 



(6) 
(7) 



For D = 0, the system decouples into two integrable 
subsystems and the dynamics is characterized by indi- 
vidual regular motions of the particle in the washboard 
potential, and bounded oscillations of the harmonic de- 
gree of freedom, respectively. For D / 0, the subsys- 
tems interact, exchanging energy. While the Q-oscillator 
performs solely bounded motion, there is the possibility 
that, for an escaping particle, the corresponding coor- 
dinate, |g| (representing the spatially-open component), 
attains large values and thus the related interaction 
forces, dW/dq and dW/dQ, vanish asymptotically, allow- 
ing transient chaos [l^-pjl]. That is, for large distance 
\q — Q\ ^> 1, the interaction vanishes asymptotically, with 
the result that the two degrees of freedom effectively de- 
couple, rendering the dynamics regular. (Although this 
decoupling is a feature of the specific example taken here, 
it is not essential for current emergence in general.) 

We focus interest on the situation when the particle is 
initially at rest at the bottom of one well of the wash- 
board potential with initial condition p(0) = q(0) = 
(representing localised initial conditions in the spatially- 
open component, situated in the centre of the interaction 
region). The total energy is initially deposited in the har- 
monic degree of freedom. The particle can then escape 
from the well only if it gains sufficient energy from its 
interaction with the excited harmonic degree of freedom. 
The total energy is fixed at E = 1.625, which exceeds 
several times the barrier height, Eb — l/ir ~ 0.318, of 
the washboard potential. A typical trajectory is shown 
in Fig. [TJ Strikingly, the particle gains sufficient energy 
from the harmonic degree of freedom that it manages to 
escape from the potential well and subsequently, avoids 
becoming trapped in wells of the washboard potential 
indefinitely. Thus, the particle moves dircctcdly to the 
right (continually increasing the coordinate, q) while the 
amplitude, Q, of the harmonic degree of freedom, per- 
forms small oscillations around the bottom of the poten- 
tial V(Q). In other words, after a chaotic transient the 
particle departs from the interaction region and the dy- 
namics of the two degrees of freedom settles on individual 
regular motion. (That the particle does, indeed, settle on 
regular motion, and that the observed behaviour is not 
merely a long-lived transient of the system, is proved in 
what follows.) 

We emphasise that the coordinate of the washboard 
particle tends to infinity and the interaction term in 
the equations of motion vanishes exponentially, exclud- 
ing the return of the trajectory to the region of the ini- 
tial conditions, as we prove below. (Remarkably, al- 
ready at t = 50 the magnitude of the interaction term 
£>tanh(<7 — Q)/ cosh(<7 — Q) has fallen below 10~ 7 .) In 
other words, the Poincare recurrence time is infinity. 
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FIG. 1: Typical time evolutions of the coordinates as indi- 
cated in the plot for D = 0.94, F = —0.5, and u = 1. Note 
the time scale on the inset figure; < t < 10 7 . 
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FIG. 2: Poincare surfaces of section represented in the (p, q)- 
plane for D = 0.94, F = -0.5, and u) = 1. 



FIG. 3: Typical time evolutions of the coordinates as indi- 
cated in the plot for D = 0.94, F = -0.5, and u = 1. The 
total energy is E = 1.625 and the boundaries of the energet- 
ically allowed region are indicated by the horizontal dashed 
lines. Superimposed is the trajectory shown in Fig. [T] 



times the period of harmonic oscillations near the bot- 
tom of a well of the washboard potential as well as of 
the potential of the harmonic oscillator. Likewise, the 
simulation time exceeds by far the time it takes for the 
particle to escape from the potential well. The Poincare 
surface of section, shown in Fig. [2j is characterized by a 
few scattered points related with the dynamics of chaotic 
transients (followed by the trajectories during the parti- 
cle's escape process) and densely covered curves associ- 
ated with the regular, rotational motion to which the 
dynamics eventually adjusts. Crucially, there are in fact 
far more trajectories evolving in the range of positive mo- 
menta than in the negative range. Hence, net motion to 
the right arises. 



In order to illustrate the dynamics of an ensem- 
ble, evolving in the four-dimensional phase space 
on the three-dimensional energy hypersurface, we in- 
voke a Poincare surface of section taken as £ = 
{q,p\Q — 0, P > }. Fig. [2] shows the Poincare surface 
of section for D = 0.94 and F = —0.5. For the computa- 
tion an ensemble of 10 5 initial momenta and coordinates 
(-P(O), Q(0)) of the harmonic degree of freedom, is taken. 
These initial values are uniformly distributed on the level 
curve E = 0.5P 2 + V(Q) + W(0, Q) in the (P, Q)-plane. 
Note that the initial conditions are unbiased, obeying the 
symmetry P <-> — P, and possess spatial symmetry with 
regard to the washboard potential. The simulation time 
interval is T = 10 5 being equivalent to almost 4 x 10 4 



It is illustrative to consider the symmetry properties 
underlying the equations of motion, Eqs. © and (J7J. 
Their structure is determined by the potential functions 
given in Eqs. (|3J)-(|5J|, which together give rise to an ef- 
fective potential U cS {q,Q) = U(q) + V(Q) + W(q,Q). 
Regarding only the washboard potential, it is periodic 
in q with period 1 and invariant under reflection in q, 
i.e., U(q) — U{—q). This amounts to symmetry with 
respect to every q n — n/2 with integer n. The inter- 
action potential is invariant with respect to changes of 
the sign of its argument, viz., (q — Q) <H> —{q — Q). In 
contrast, the tilted harmonic oscillator potential V(Q) is 
without reflection symmetry and, consequently, the sys- 
tem of Eqs. ([5]) and (J?} is not invariant under inversions 
Q ■(->• — Q. The symmetry lines are given by 
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Si 

s 2 



cosh((7 — Q) 

Q + F-^ ta y-^ =o. 

cosh(g — Q) 



(8) 
(9) 



The symmetry line S\ exhibits the following symmetry: 
Q-+-Q, -+q^---q: Si >n -+-Si- n 

( 10 ) 

with integer n labelling the several branches of the sym- 
metry line as Si )7l . A single branch, S 2 , arises from the 
second equation with no apparent symmetry as shown in 
coordinate space in Fig. [3] The branches of the symme- 
try lines Si, n related to the starting potential well and 
its neighbour to either side are displayed in the figure, 
i.e., those with n = —2,-1,0,1,2. Crucially, reflections 
in the (q, Q)-coordinate plane on the symmetry line S 2 
do not respect inversion symmetry, and thus, no pairs of 
counterpropagating trajectories are produced by the lo- 
calised, but otherwise unbiased, set of initial conditions. 
Symmetry breaking of the trajectory displayed in Fig. [T] 
becomes apparent by superimposing the latter on the di- 
agram in coordinate space in Fig. |3J In particular after 
an ultimate crossing of the symmetry line Si t2 the trajec- 
tory proceeds exclusively to the right as indicated by the 
arrow. This spatial symmetry breaking is a prerequisite 
for the occurrence of directed net motion. In conclu- 
sion, even though the initial conditions are unbiased for 
the ensuing dynamics, the time- reversibility symmetry is 
broken. 

We now prove that there exist trajectories which, once 
they have entered a certain region in phase space, provide 
continual transport of washboard particles. The Hamil- 
tonian associated with the system in Eqs. ©,([7]) reads 
as 



H 



1 , „, , 1 



2 1 TTI„\ I 



U(q) 



V(Q) + W(q,Q) 



T w 2 
= H 1 (p,q) + H 2 (P,Q)+W(q,Q), 



(11) 



with U, V, W given in © , g]) , © . 

Theorem: Consider the Cauchy problem 
for the Hamiltonian (ITTI) with initial data 
(p(0),<z(0),P(0),Q(0)). 

Let E denote the total energy of the system and let Q 
denote the positive root of the equation V(Q) = E, i.e., 
the maximal coordinate that the oscillator could attain if 
all of the energy were in the oscillator degree of freedom. 

Assume that the energies of the washboard particle 
and of the harmonic oscillator satisfy the constraints 



Ei(0) >E a + D 



1 



cosh(<7(0) - Q) 



and 



£ 2 (0) < E — D — E a , 
respectively, where 

E a >-, 



(12) 
(13) 
(14) 



the latter being the energy barrier height of the wash- 
board degree of freedom. 

If, in addition, the intial values of the coordinates of 
the washboard and oscillator satisfy the constraint 

?(0)-Q(0) >g(0)-Q>sinh- 1 (l), (15) 

then it holds that, for some p, 

p(t)>p>0, (16) 

for all times t > 0. 

Proof: The proof utilises that existence of solutions 
with p > is guaranteed by the asymptotic behaviour of 
the washboard coordinate as q —> 00. The momentum p 
is represented as follows: 



1 



p 2 =E-U(q)-E 2 -W(q,Q). 



(17) 



With the help of (Til?)) we obtain the lower bound for the 
initial value 



1 



p 2 (0) = E a -max.[U(q(p))]+D 



1 



1 1 

> °--+ £) cosh ( g ( )_Q( )) 



cosh(g(0) - Q(0)) 
(18) 



From the lower inequality (fT4"l) follows that the momen- 
tum p(0) is positive, regardless of the value of the positive 
term D/ cosh(q(0) - Q(0)). 

The change in energy of the two subsystems under the 
dynamics is determined by 

§ = « 

Suppose (in order to reach a contradiction) that the tra- 
jectory reaches a point at which p — 0, i.e., that the 
trajectory reaches a stationary point in the washboard 
coordinate, q, after which it might turn back. Let q\ > qo 
denote the first point at which that happens. Using that 
p = dq/dt and P — dQ/dt one obtains 



w m jp n f 1 tanh(g - Q) 

m)-E 2{ o)=D[ Ql ^^- dQ . 

Jqo cosh(g-Q) 



(21) 
(22) 



The maxima of the integrand tanh(g — Q)/cosh(q — Q) 
are located at \q — Q\ = sinh _1 (l) and for 



Iff — Q| >sinh- 1 (l) 



(23) 



the integrand remains non-negative. Hence, provided the 
coordinates stay in the range given by (|15j) the washboard 
degree of freedom loses energy while the harmonic oscil- 
lator gains energy. We derive an upper bound on the 
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magnitude of the energy loss of the washboard particle 
and demonstrate that, provided the above assumptions 
on the initial data hold, this loss cannot in fact be suffi- 
cient to halt or turn back the particle. To this end, the 
integration in (|2 1[) is performed yielding 



#!(!)- #!(()) 



-d r tanh f g -Qu 



> -D 



= -D 



> -D 



qo cosh(q-Q) 
qi tanh(g - Q) 
qo cosh(q - Q) 
1 



(24) 



dq 



1 



cosh(gi — Q) cosh(<?o — Q) 
1 



cosh(q - Q) 



(25) 



Suppose that initial data for the coordinate q(0) satisfies 
the inequality (fT5"j) and, in addition, the washboard de- 
gree of freedom satisfies our energy constraints (fT2"|) and 
(H2), then 



Ei(l) > Ei (0) — D 
> E a >-. 



1 



cosh(q - Q) 



(26) 
(27) 



Thus, the energy of the particle in the washboard degree 
of freedom remains above the separatrix level, contra- 
dicting the supposition that the particle halts (attaining 
p = at some position qi), and we therefore have that 
p(t) > for all times t > 0. Indeed, since the energy 
E\ of the washboard particle is bounded away from the 
barrier level, 1/vr, it follows that p(t) > p > for some 
p, which implies that q — !> oo as t — !> oo. □ 
That regions of phase space satisfying the assumptions 
of the theorem do, indeed, exist, will be shown below. 
Firstly, we consider the asymptotic behaviour of such a 
transporting trajectory. Given that q — > oo, the asymp- 
totic energy of the washboard particle is bounded by 



£i(oo) > Ei(0) 



D 



cosh(<jo - Q) 



(28) 



With the assumption in Eqs. (|I2I) and (IT41 we deduce 
that 



Ei{oo) >E a >- 

7T 



(29) 



The last property means that the washboard particle 
loses only a finite amount of energy and the energy of 
the washboard particle approaches the level E = E a from 
above. At the same time the energy of the harmonic os- 
cillator asymptotically approaches E2 = E — E a — D 
from below. Hence, the coordinate of the harmonic os- 
cillator performs bounded oscillations around Q = 
which, together with growing q, assures that the distance 
q(t) — Q(t) effectively grows as time progresses. Hence, 




60 



80 



100 



FIG. 4: Time evolution of the energy of the washboard parti- 
cle Ei and the harmonic oscillator E%. The two dashed hori- 
zontal lines correspond to the energy levels 1/tt and E—D—E a 
bounding the energies Ei = E a and Ei — E — D — E a in the 
asymptotic regime. 



for initial coordinates fulfilling the constraints of the the- 
orem, the integrand in (|24p vanishes in an exponential 
fashion. 



In conclusion, as soon as the inequalities (fl2^ .([T3 j) 
and (|I5[) . imposing conditions in phase space, are satis- 
fied, the process of redistributing the finite (albeit small) 
amount of energy DJ cosh(q — Q), that the washboard 
degree of freedom possesses in excess of the asymptotic 
level E a , into the harmonic oscillator is irreversible. This 
asymptotically decouples the two subsystems and, most 
importantly, the washboard particle, being equipped 
with energy Ei > E a > 1/tt, moves in an unhindered 
and unidirectional manner, which establishes continual 
transport that cannot be terminated. In Fig. 21 we illus- 
trate a trajectory that shows an initial energy redistri- 
bution from the harmonic oscillator into the washboard 
particle, after which the phase space region conforming 
to the assumptions of the theorem is reached. Notice the 
pronounced growth of the energy Ei at the expense of 
Ei for t > 25. Crucially the energy E2 falls below the 
level E — D — E a at t ~ 27, that is at the moment of 
no return. Strikingly, this happens at the moment when 
the difference q — Q starts to increase beyond the value 
sinh _1 (I). Afterwards the above mentioned (small) en- 
ergy loss (gain) in Ei (E2) takes place as a result of which 
the energy Ei settles onto E a . The latter value lies above 
the lower bound I /tt as given in (|I4|) (represented by the 
lower dashed line in Fig. |4|) while the energy E2 attains 
an asymptotic level below E — D — E a (represented by 
the upper dashed line in Fig. [4]). 

In summary, we have demonstrated that it is possible 
to obtain directed motion in autonomous Hamiltonian 
systems under the minimal conditions that (i) transient 
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chaos is supported and (ii) some of the degrees of free- 
dom serve to break the time reversibility with regard to 
a set of unbiased localised initial conditions. (The model 
system is minimal since these two conditions constitute 
the only indispensable prerequisites. Notice in particular 
that no additional modulation held is required in order to 
break the necessary spatio-temporal symmetries.) That 
chaos is needed only in an initial stage of the dynamics 
in order to guide the trajectory into the range of regular 
motion has a drastic implication, namely the directed net 
motion is provided by regular motion and thus is of con- 
tinual nature. This is at variance with previous studies 
of transport in non-autonomous [3 - [l] and autonomous 
Hamiltonian systems [s| , [nj , where a mixing phase space 
is needed to support directed chaotic transport. Despite 
the fact that finite asymptotic currents were observed, 



sustaining transport in general in nonintegrable systems 
seems to be hampered because of the intermittent char- 
acter of the dynamics inducing power-law decay statis- 
tics of the duration of periods of regular motion. On 
the other hand, provided the minimal prerequisites men- 
tioned above are given, application of our approach to 
accomplish directed net motion in other Hamiltonian sys- 
tems, such as infinite lattice systems, is straightforward. 
For example, in the context of the Holstein model, the 
charge being initially trapped in a confined region of the 
molecular chain constitutes localised initial conditions. 
Charge motion along a molecular chain can be directed 
by taking into account non-reversion symmetric oscilla- 
tors representing the intra-molecular vibrational degrees 
of freedom. 
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